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ABSTRACT 

We use absolutely calibrated data between 3 and 90 GHz from the 2006 balloon flight of the AR- 
CADE 2 instrument, along with previous measurements at other frequencies, to constrain models of 
extragalactic emission. Such emission is a combination of the Cosmic Microwave Background (CMB) 
monopole, Galactic foreground emission, the integrated contribution of radio emission from external 
galaxies, any spectral distortions present in the CMB, and any other extragalactic source. After re- 
moval of estimates of foreground emission from our own Galaxy, and the estimated contribution of 
external galaxies, we present fits to a combination of the flat-spectrum CMB and potential spectral 
distortions in the CMB. We find 2 a upper limits to CMB spectral distortions of fi < 5.8 x 10 -5 and 
| Iff | < 6.2 x 10 -5 . We also find a significant detection of a residual signal beyond that which can be 
explained by the CMB plus the integrated radio emission from galaxies estimated from existing sur- 
veys. After subtraction of an estimate of the contribution of discrete radio sources, this unexplained 
signal is consistent with extragalactic emission in the form of a power law with amplitude 1.06±0.11 K 
at 1 GHz and a spectral index of —2.56 ± 0.04. 

Subject headings: cosmology: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) is cur- 
rently our most precise window on the physics of the early 
universe. Measurements of the frequency spectrum of the 
CMB can rule out alternative cosmologies and place lim- 
its on physical processes that may distort the spectrum, 
including dark matter particle decay and reionization. 
Departures from a thermal blackbody spectrum are ex- 
pected at a small level from a variety of mechanisms. 

The Cosmic Background Explorer (COBE) satellite 
observed the spectrum of the CMB with the Far- 
Infrar ed Absolute Spectr ophotometer (FIRAS) instru- 
ment (|Mather et al.|[l990l ) at wavelengths between 1 cm 
and 1 00 jim. FIRAS re sults reporte d by iFixsen et alJ 
1996ft . iMather et all (| 19991 ) and IFixsen fc Mather! 



2002) are consistent with a blackbody spectrum at a 
temperature of T C mb = 2.725 ± 0.001 K. 

Absolutely calibrated measurements of the CMB 
at longer wavelengths (lower frequency) than FIRAS 
have been performed with ground-based and balloon- 
borne experiments. Among the most sensitive of these 
measu re ments are t hose of _ Johnson and Wil kinsonI 
(11987ft. iLevin et alJ 



Bersanelli et al 
Stages et al 



(11992ft . iBersanelli et all (11994ft . 
-19951), Istaggs et all (l!996al). 
"(Il996bl) , iRaghunathan fc Subrahmnavanl 



btaggs et al. (lyyobj, Kagnunatnan & bubranmnavan 
(j2000ft . IFixsen et alJ (|2004h . iSingal et alJ (|200l . and 
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IZannoni et ail (|2008ft . 

The second generation of the Absolute Radiometer 
for Cosmology, Astrophysics, and Diffuse Emission (AR- 
CADE 2) was conceived as a balloon-borne experiment to 
improve constraints on spectral distortions in the CMB, 
with particular emphasis on the 3 — 10 GHz frequency 
range. ARCADE 2 uses a unique, clear aperture instru- 
ment design with the bulk of the instrument operating at 
or near the temperature of the CMB. This minimizes the 
potential contribution to instrument systematics from 
warm, emissive optics. The instrument uses a set of mi- 
crowave feed horns to compare the sky to a large, cryo- 
genic blackbody calibration target. The results described 
here are from the se cond version of the instrument, de- 
scribed in detail by ISingal et al.l (|2008ft . The sky mea- 
surements fr om the secon d flight of this instrument are 
described bv lFixsen et al.l (|2008j ). which presents a detec- 
tion of extragalactic emission consistent with a power law 
plus constant CMB temperature. The model of Galac- 
tic emission used in interpreting the ARCADE 2 data is 
described bv lKogut et"ai1 (|2008ft . 

In this paper, we use the combination of ARCADE 2 
and other data sets to present a detection of 3 GHz 
emission in excess of that expected from the CMB 
and existing source counts of radio galaxies. The ex- 
cess emission amplitude and spectral index described 
here differ fr o m the extragalactic emission described by 
IFixsen et all (|2008ft . Here, we are concerned with the 
residual excess emission that cannot be explained from 
known classes of emission, and explicitly correct for the 
estimated contribution of discrete radio galaxies. We also 
use the combined data to place limits on spectral dis- 
tortions to the CMB, and show that canonical spectral 
distortions cannot explain the excess emission. 

This paper is organized as follows: Section 2 summa- 
rizes the estimates of isotropic, extragalactic emission at 
a variety of frequencies that we have used in our analysis. 
Section 3 examines the potential contribution of extra- 
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galactic point sources and their potential to affect our 
conclusions. Section 4 presents our spectral fits to the 
data and our limits on spectral distortions of the CMB. 
Section 5 presents discussion of the results, including po- 
tential explanations for the source of the excess emission. 

2. RESULTS FROM ARCADE 2 AND OTHER 
SURVEYS 

For our analysis, we use the data from the 2006 flight 
of the ARCADE 2 instrument, from FIRAS, and from 
lower frequency ground-based surveys. FIRAS measures 
a high-precision difference between the sky and a cali- 
brated reference target. The result is a set of values with 
tiny relative errors, and a larger, 1 mK calibration error 
common to all the data points. Table [T] summarizes the 
remaining data used in our anal ysis, which include s AR- 
CADE 2, the 22 MHz survey of iRoger et alJ (| 1999( 1. the 
45 M Hz survey oflMaeda et alJ (|1999f) . the 408 MHz sur- 
vey of | Haslam et al.l (Il981|). and the 1.42 GHz survey of 
iReich k Reichl (|1986pKogut et alJ (|2008h describes the 
process of estimating the Galactic component from each 
of these data sets. The data in TableQ]is the resulting es- 
timate of the residual, extragalactic, isotropic emission. 
The ARCADE 2 data in the 3 - 10 GHz range are shown 
in Figure 1; they lie significantly above the 2.725 K black- 
body CMB determined by FIRAS at higher frequencies. 
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Fig. 1. — Detection of extragalactic radio emission by AR- 
CADE 2 beyond the contribution of discrete radio sources 
and the expectation of 2.725 K blackb ody radiation . Dat a 
points are the ARCADE 2 results from Fix sen et alJ (120081 ). 
and have be en corrected for M ilkv Way Galactic emission 
described bv lKogut et alJ (|2008l ). The dashed curve is a con- 
stant 2.725 K blackbody, consistent with FIRAS measure- 
ments of the CMB. The dot dash curv e is an estimate of the 
discrete radio source contribution from Gcrvasi et al. (2008a) 
model "Fitl" added to the 2.725 blackbody. The data points 
lie significantly above this dot dash curve, indicating our de- 
tection of unexplained, excess emission. The solid curve is 
the best fit of the combined data of Table 1 and FIRAS to a 
power law plus a constant CMB temperature. 



In ou r analysis, we have e xclud ed the 100-200 MHz re- 
sults of |Rogers&iBowma^ (JUS! . They find a minimum 
diffuse background of 237 K at 150 MHz, but their work 
does not provide an independent estimate of the Galac- 
tic contribution. We can, however, check for consistency 



by us ing the Galactic model described by iKogut et al.1 
(2008) extrapolated to 150 MHz, where we find an ap- 
proximately 60 K Galactic contribution to the diffuse 
background in the region of high Galactic latitude. Ap- 
plying this correction, we find both the emission ampli- 
tude and spectral index are consistent with the fits we 
present in section 4. 

We have not included a number of other mea- 
su rements, including the roc k et-bor ne measurements 
of iGush. Halpern. fc Wishnowl (|1990f) and the ground- 
based and balloon-borne measurements cited earlier. 
The size of the uncertainties quoted in these measure- 
ments results in no significant contribution to the con- 
straints on our model fits. 

3. CONTRIBUTION OF SOURCES 

The set of measurements in Table 1 do not have suffi- 
cient angular resolution to reject discrete radio sources. 
Instead, we must estimate the contribution of these 
sources through one of two ways: direct radio surveys 
designed to detect such sources, or measurements of the 
far-IR background which can trace the integrated emis- 
sion of such sources through the correlation of the far-IR 
and radio emission. We examine these two methods in 
turn. 

3.1. Expectation from source counts 

The sky brightness temperature contributed by dis- 
crete sources can be composed as the sum of two parts: 
the source population that has been characterized by ex- 
isting surveys and the contribution of sources below the 
flux limit of existing surveys. We write this as: 



T = T(S > S Umit ) + 



2fc R 



'limit dN 

~d~S 



SdS, (1) 



where T(S > Slimit) is contribution from sources with a 
flux S greater than the survey limit Slimit- The wave- 
length of observation is A and k-Q is the Boltzmann con- 
stant. We characterize sources below the survey limit 
with their differential number counts, dN/dS, and as- 
sume that there is a lower limit cutoff to the source pop- 
ulation at a flux of S mm . Radio source count surveys 
reveal a faint source population with differential number 
counts proportional to a power law 



dN 
~d~S 



= S~ 



(2) 



An index 7 of 2.5 corresponds to a static, Euclidean uni- 
verse with uniform filling of sources, whereas faint radio 
surveys find 7 in the range of 2.0 to 2.3. Such source 
counts can not extend to arbitrarily low fluxes, or the 
total contribution would diverge. A realistic distribution 
of sources, of course, would not have a sharp cutoff at 
S m i„. In practice, we can characterize S m i n as the flux 
below which the index 7 falls below 2, as there will be 
negligible additional contribution to the integral below 
this limit. 

Deep surveys of radio sources have been performed 
at a number of frequencies. Particularly useful are the 
surveys at 1.4 and 8.4 GHz w ith the Very Large Array 
(VLA). iFomalont et alJ (|2002t ) reports the results of an 
8.4 GHz survey with a survey limit of 7.5/zJy and finds a 
faint-end index to the number counts of 7 = 2.11 ± 0.13. 
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TABLE 1 

Measurements of Extragalactic Radio Emission 



Frequency (GHz) 


Extragalactic 3 Temperature (K) 


Error (K) 


Radio Sources' 3 


Error c 


Residual Emission" 1 


Error 


0.022 


21200* 


5125 f 


7090 


709 


14100 


5175 


0.045 


4355 g 


520 s 


1020 


102 


3334 


530 


0.408 


16.24 h 


3.40 h 


2.61 


0.26 


13.6 


3.41 


1.42 


3.213 j 


0.53 j 


0.089 


0.009 


3.124 


0.53 


3.20 


2.792 k 


0.010 k 


0.010 


0.001 


2.706 


0.010 


3.41 


2.771 k 


0.009 k 


0.008 


< 0.001 


2.682 


0.009 


7.97 


2.765 k 


0.014 k 


< 0.001 


< 0.001 


2.577 


0.014 


8.33 


2.741 k 


0.016 k 


< 0.001 


< 0.001 


2.545 


0.016 


9.72 


2.732 k 


0.006 k 


< 0.001 


< 0.001 


2.505 


0.006 


10.49 


2.732 k 


0.006 k 


< 0.001 


< 0.001 


2.488 


0.006 


29.5 


2.529 k 


0.155 k 


< 0.001 


< 0.001 


1.887 


0.151 


31 


2.573 k 


0.076 k 


< 0.001 


< 0.001 


1.900 


0.074 


90 


2.706 k 


0.019 k 


< 0.001 


< 0.001 


1.098 


0.015 



a This is the monopole temperature with the Milky Way Galactic contribution removed as by[Kogut ct al] (120081 ) 

b Estimate of extragalactic discrete radio source contribution from Gcrvasi ct al. (2008a) model "Fitl". Units are K, antenna 

temperature. 

c We have adopted a 10% fractional error for the[G ervasi et al.l (12008a!) "Fitl" model (see text). 

d Residual extragalactic emission after subtraction of radio source. Values have been converted to K, antenna temperature. 
c Error in residual extragalactic emission after subtraction of radio source, K, antenna temperature. 

f Data from [Roger et aTl ||1999D corrected for Galactic emission with the model described by IKogut et al.l d2008l) . Units are K, 

antenna tem perature. 

g Data from Maeda ct al. (1999) corrected for Galactic emission with the model described by IKogut et al.l ( 1200811 . Units are K, 
antenna temperature. 

h Data from Haslam ct al. (1981) corrected for Galactic emission with the model described by[Kogut ct al. (2008]). Units are K, 
antenna tem perature. 

j Data from Reich & Reich (1986) corrected for Galactic emission with the model described by[Kogut et al. (2008f). Units are K, 
antenna temperature. 

k ARCADE 2 IFixsen etaLl (120081 ). Units are K, physical temperature. 



iWindhorst et~ai] (l993) argue that the sources in the nJy 
flux range are dominated by ordinary spiral galaxies, 
which produces a natural lower limit of 30 nJy; below 
this limit there are insufficient galaxies. 

Figure 2 shows a range of estimates for the contribu- 
tion of discrete sources to the ARCADE 2 3.20 GHz mea- 
surement. We have calculated the expected contribution 
to the extragalactic sk y temperature using the results of 
iFomalont et alj ((2002) and varying the faint end index, 
and plotting as a function of S m m- We have scaled the 
temperature from 8.4 to 3.2 GHz using a frequency spec- 
tral index of -2.75. This spectral index is characteristic 
of starburst and normal spiral galaxies with synchrotron 
emission, though at 8.4 GHz there could be some con- 
tribution from flat spectrum sources suc h as AGN. The 
source s fainter than 35 /zJy described bv lFomalont et aTl 
(2002) have a spectral index distribution that peaks at 
-2.75. From this analysis, we conclude that the contri- 
bution is likely in the range of 5 to 10 mK at 3.20 GHz. 

It is interesting to consider to what extremes one would 
need to take the number counts in this analysis to ac- 
count for the measured 58 mK excess at 3.20 GHz. Using 
an in dex of 7 = 2.5 and ext rapolating the number counts 
from [Fomalon t et al.l (|2002D , we find that we would need 
to extend the lower flux limit to 0.3 nJy, which would 
result in a source density of ~ 8 x 10 per square arcmin. 
This extreme scenario is not a plausible distribution of 
ordinary galaxies. We also note that the normalization 
of the differential number counts in the 8.4 GHz survey is 
sufficiently accurate to not contribute a significant source 
of error to this analysis. 

We have focused on the 8.4 GHz survey results, as these 
measurements are amenable to extrapolation to fainter 
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Fig. 2. — Estimated contribution of extragalactic source 
counts to the sky brightness at 3.20 GHz, versus the assumed 
faint end cutoff of 8.4 GHz source counts. The sky brightness 
is first calculated at 8.4 GHz from the sum of the existing 
source population greater than 7.5 ^iJy (Fomalont et al 2002) 
and the contribution of fainter sources assuming a continua- 
tion of the differential number counts with a power law index 
of 2.1 (solid curve), 2.2 (dotted curve), 2.3 (dashed curve), 
2.4 (dash dot curve), and 2.5 (dash dot dot curve). The con- 
tribution is scaled to the ARCADE 3.20 GHz channel with a 
frequency index of -2.75, as is typical of sources in faint radio 
surveys. The measured index of differential number counts 
at the faint end of the 8.4 GHz survey is 7 = 2.11 ± 0.13. 

fluxes and to the ARCADE 2 frequencies. Our results, 
however, do not d epend on the results fro m a radio survey 
at one frequency. iGervasi et al.l (|2008a|) describe a more 
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comprehensive analysis of the potential contribution of 
unresolved extragalactic radio sources by examining data 
from a wide variety of radio surveys, and fitting an em- 
pirical, two-population model to the survey data at each 
frequency. Their combined result for extragalactic radio 
source brightness versus frequency is described by a sin- 
gle power law model, "Fitl", with amplitude 0.88 K at 
0.61 GHz, and a power law index of -2.707. This model is 
also shown in Figure 1, indicating again that this contri- 
bution is insufficient to explain the ARCADE 2 results. 
The brightness temperature values of this model, evalu- 
ated at the frequencies used in our analysis, appear in 
Table 1. Because the sample size and brightness lim- 
its of the radio surveys vary with frequency, the error 
in their estimate also varies with frequency. These fre- 
quencies are not identical to the frequencies measured by 
ARCADE 2. For Table 1, we adopt a fractional error es- 
timate of 10% independent of frequency, w hich we believe 
is a con servative estimate of the error in the lGervasi et all 
(2008a) modeling. These estimates for the contribution 
of extragalactic sources are consistent with the analysis 
described above. 

It also is interesting to consider whether the existing 
surveys have missed significant flux from the sources. 
The low frequency faint source observations are primar- 
ily interferometric and have the possibility of overresolv- 
ing the source and missing flux in extended lo w surface 
brightness emission. iHenkel fe Partridge! (|2005h consider 
the evidence for this and conclude that 20% may be an 
upper limit to t h is effe ct for mJy flux levels at 8.5 GHz. 
iFomalont et alj (|2006| ) suggests that at 1.4 GHz only 
a few percent of sources are larger than 4 arc seconds 
and that other reports of a larger figure in other sur- 
veys are actually con fusion of multiple disparate sources. 
iGarrett et all (2000) compare their 1.4 GHz survey con- 
ducted using the Westerbork Synthesis Radio Telescope 
and its larger, 15 arcsec effective beam with previous 
VLA measurements of the same region with a similar 
noise level. Of a total of 85 sources in their survey, they 
find 22 not apparent in the previous VLA survey. Some 
of these 22 likely correspond to the combined flux of mul- 
tiple sources that were resolved by the VLA measure- 
ments. At least 2 sources, however, appear to be rela- 
tively nearby discrete sources with emission from a large 
enough region to have been resolved out by the VLA mea- 
surements. We conclude from these studies that it seems 
unlikely that sufficient flux has been missed in surveys of 
known objects to explain our residual emission. 

Another method to examine the possibility of ex- 
tended low surface brightness emission in extragalac- 
tic sources is ra dio observations of the halo s of n earby 
edge on spirals. Ilrwin. English, fe Sorathial (|1999D and 
Ilrwin. Saikia. fc Englisbl (|2000D report results of VLA 
surveys for radio emission from nearby edge-on spi- 
rals. These studies can elucidate the connection be- 
tween the star formation processes that drive the far- 
IR background, and the supernova processes that drive 
radio emission, but do not provide evidence that large 
amounts of radio flux are missed in surveys of more dis- 
tant sources. 

3.2. Connection with Far-IR Background 

The cosmic far-IR background has been detected at 
a level of approximately 10 — 20nWm~ 2 sr~ 1 with FI- 



RAS and DIRBE flPuget et alj |1996t iFixsen et al.1ll998l ; 
lHauser et al.l Il998h . We can use the univers al radio 
to fa r-IR correlation in star-forming galaxies (|Condonl 
1992) to estimate the expected extragalactic radio back- 
ground that can be attributed to galaxies contributing 
to the cosm ic far-IR background. [Haarsma fc Partridge! 
(|2000f) and iDwek fc Barkerl (|2002j) specifically address 
this prospect. The conclusion of these studies is that 
the far-IR measurements are consistent with the existing 
surveys of radio galaxies described earlier. For example, 
the radio brightness temperatur e of 18 K at 178 MHz pre- 
dicted bv lDwek fc Barkerl (|2002h i s within 1 a of t he radi o 
galaxy contribution modeled by iGervasi et al.l ((2008a) . 
This emission is insufficient to account for the excess de- 
tected by ARCADE 2. 

There is no obvious way around this limit by consider- 
ing departures from the far-I R radio correlation associ - 
ated with faintness or redshift. iGarn fc Alexander! (2008) 
stack IR-selected galaxies and data from faint radio sur- 
veys and find there is no evidence for a change in the 
far-IR to radio correlation with fainter galaxies. Simi- 
larly, there does not appear to be evidence fo r a change 
in the far-IR to radio cor relation with redshift (jChapmanl 
l2005t lFraver"et~aT1[2006Tl . 

There is some room for breaking the far-IR to ra- 
dio correlation; the physical processes of far-IR emission 
from dust heated by star formation and radio emission 
driven by supernovae are relate d but differ in timescale 
(see, e.g., I Murphy et al.l 120081 ). There does not ap- 
pear, however, to be a case for a sufficient difference 
in timescale to account for our measurements in known 
populations of galaxies. 

4. FIT OF MODEL SPECTRUM 

Here we fit a model spectrum to the combined data 
set of the residual emission column of Table 1 and the 
higher frequency results from FIRAS. The form of the 
fitting function is: 

T{y) =T +A(v/1 GUzf + AT(i/), (3) 

where To is the CMB baseline temperature, A is the 
power law amplitude at 1 GHz, f3 is the power law index, 
v is frequency, and AT(v) is a CMB spectral distortion. 
Data are converted (where necessary) to units of antenna 
temperature before fitting using: 

hv jk 



u/llz/fcTphys 



r 



(4) 



where h is Planck's constant, k is Boltzmann's constant, 
T Aat is antenna temperature, and T PllyB is the physi- 
cal temperature. We use a Levenberg-Marquardt non- 
linear lejst^c£u^e^_minimization for the fitting proce- 
dure (Marquardt 1963). We have allowed for the 1 mK 
relative calibration uncertainty between the FIRAS tem- 
perature scale and the other measurements in determin- 
ing the size of the errors on the fit parameters. This 
was accomplished by adding or subtracting 1 mK to the 
FIRAS measurements and noting the change in the pa- 
rameter value. This change was then added in quadra- 
ture to the error derived without allowing the relative 
calibration error. In most cases this makes a only very 
small difference. The exception is Tomb, where we re- 
cover the expected 1 mK error. The parameters for the 
fits described in this section are summarized in Table 2. 
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4.1. Power-Law Plus CMB 

Figure 1 shows that there is clear excess emission de- 
tected by ARCADE 2 in the 3 GHz channels compared 
to what is expected from the CMB plus the contribution 
of extragalactic radio sources. The unexplained resid- 
ual emission is consistent with a power law with ampli- 
tude 1.06 ±0.11 K at 1 GHz, with a spectral index of 
)3 = -2.56 ± 0.04. 

We have also experimented with inflating the assumed 
errors for the removal of extragalactic discrete sources. 
As noted in Section 3.1, we have assumed a fractional 
error of 10%, independent of frequency, for the discrete 
source contribution model. Inflating this error to 50% 
fractional error makes much less than a 1 <r change in 
the values of our power law fit parameters above, and 
only a small (~ 10%) increase in the quoted errors in the 
fit parameters. This is because the error in the removal 
of the discrete sources is smaller than the other errors 
in the low frequency measurements, as can be seen by 
inspection of Table 1. 

4.2. Free- Free Distortions 

Free-free distortions to the CMB spectrum 
can arise from energy released at lower redshifts 
(IBartlett fc Stebbinsl[l99lt iGnedin & Ostrikenll997l : lObl 
1999), and can be characterized by 

ATg(u)=T ^, (5) 

where Yg is the optical depth to free-free emission, T is 
the undistorted CMB temperature, x is the dimension- 
less frequency hv/kTo, and AT is apparent temperature 
distortion. 

A lower limit can be placed on the optical depth to 
free-free emission from late time effects ofYff>8xlO~ 8 
(|Haiman fc Loe b 1997). The current upper limit of 
Yg < 1.9 x 10~ 5 comes from combining data from FI- 
RAS and previous ground-based CMB spectrum mea- 
surements (|Bersanelli et al.lll994ft . 

We have performed a four component fit to the data 
to assess whether there is evidence to support a free-free 
spectral distortion to the CMB, compared to the three 
parameter fit described in the previous section. The four 
fit components consist of a constant CMB temperature, 
a power law amplitude, a power law index, and a free-free 
amplitude. The fit parameters are presented in Table 2. 
The addition of the free-free amplitude does not improve 
the reduced \ 2 OI the fit compared to the fit presented 
in Figure 1 and is therefore not justified by the data. 

We have also examined a two parameter fit consisting 
of a constant CMB component and free-free distortion 
component; the parameters are shown in Tabic 2. This fit 
produces a significantly worse reduced % 2 , and is there- 
fore not consistent with the source of the unexplained 
emission. 

The 2a limits on the free-free amplitude at 1 GHz de- 
rived from our four parameter fit are — 0.44K < ATy {t < 
0.52K. This corresponds to an upper limit on the free- 
free optical depth of 

lYffl < 6.2 x 10" 5 (6) 

This limit is le s s con strai ning than those repor ted by 
iBersanelli etail (|1994l) and lGervasi et all (|2008bD . This 



is the result of the additional degrees of freedom al- 
lowed by our four parameter fit. As described above, 
the four parameter fit is a better description of the data 
than the CMB p l us fre e- free distortion fits perform ed by 
IBersanelli etal) (|1994D and iGervasi et all (|2008bP . We 
conclude that the tighter constraints offered by those 
studies are likely too optimistic. 

It is interesting to consider how future measurements 
might improve the Yg limit. We have used the exist- 
ing fits and asked how much tighter the limit becomes if 
an additional measurement of some fixed fractional ac- 
curacy is added to the data set. We have run this test 
as a function of frequency; the frequency range with the 
greatest effect is between 0.3 and 3.0 GHz. where a fac- 
tor of several tighter constraint is potentially achievable. 
The results are shown in Figure [3j where we have plot- 
ted the size of the 2 a errors on If as a function of the 
frequency of the additional measurement. We have as- 
sumed that the measurement point falls precisely on the 
existing fit. Better fractional error results in tighter con- 
straints on Yg. We have also considered the impact of 
repeating this test with additional measurements at 3, 5, 
8, 10, 30, and 90 GHz with 0.002 K precision, as might 
be obtained with a future flight of ARCADE. Such a 
measurement would enhance the power of future lower 
frequency measurements to constrain Iff. 
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Fig. 3. — Improvement of limit on Yg obtainable with future 
measurements. Shown are the sizes of the 2 a errors on Yg 
(red curves) that would result from an additional measure- 
ment at a frequency shown on the x-axis, with a precision of 
0.1% (dotted curve), 1% (dashed curve), or 10% (dot dashed 
curve). The blue curves repeat this test under the assump- 
tion that there are additional measurements at 3, 5, 8, 10, 30, 
and 90 GHz with 0.002 K precision, as would result from a 
future flight of ARCADE 2. The additional measurement is 
assumed to be centered on the value of the existing fit. 

4.3. Chemical Potential (fj,) Distortions 

As we have discussed, energy release early in the 
universe can distort the spectrum of the CMB. En- 
ergy injection after a redshift of ~ 10 6 no longer re- 
sults in a planckian spectrum, but instead forms a spec- 
trum with a Bose-Einstein photon occupation number 
dSunvaev fc Zeldovichlll970l) : 

= e s+M(*) - 1 ' (7) 
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TABLE 2 

Spectral Fits to Combined ARCADE 2, 



FIRAS, Radio Survey Data 



Parameter 


Power Law Fit 


Power Law + Yg 


Y ff Only 


Power Law + fi Distortion 


TCMB a 


2.725 ± 0.001 


2.725 ± 0.001 


2.725 ±0.001 


2.725 ±0.001 


Power Law Amplitude b 


1.06 ± 0.11 


1.00 ±0.37 




1.05 ±0.11 


Power Law Index 


-2.56 ±0.04 


-2.58 ±0.11 




-2.56 ± .05 


Free-Free Amplitude c 




0.04 ± 0.24 


0.54 ±0.07 




fj, Amplitude 








(-0.73 ± 3.3) x 10- 5 


Degrees of Freedom 


53 


52 


54 


52 


x 2 


60.8 


60.8 


107.1 


60.8 


Reduced \ 2 


1.15 


1.17 


1.98 


1.17 



a The best fit thermodynamic temperature of the CMB in K. 

b The fit amplitude of a power law component in K (antenna temperature) at 1 GHz. 
c The fit amplitude of ATy ff in K (antenna temperature) at 1 GHz. 



where x = hu/kT is the dimensionless frequency, 
and fJ-(x) is a frequency-depend ent chemical po- 
tentia l. A series of papers fDanese fc De Zotta 
1980t fBurigana. Danese. & De Zottil Il991allbl ; 
Burigana. De Zotti. k Danesd I1995T T has investigated 
in detail the shape of the resulting spectral distortions 
after inclusion of free-free processes which act at the low 
frequency end. The shape of the distortion depends on 
the range of redshift over which such energy injection 
takes place. We have used their analytic description of 
such distortions to the CMB to provide a functional 
form to fit the data in Table 1. A necessary ingredient 
is the baryon density, for which we adopt a value of 
Q b h 2 = 0.02265, from the 5-year WMAP data, includ- 
ing constraints from Baryon Acoustic Oscillation and 
Super nova measurements as described by Hinshaw et al. 
((20081 ). 

In addition to the [i distortion, we have included a 
power law amplitude and index to the fit parameters. 
Figure 4 shows the result of the fit as well as the residuals. 
The addition of the [i distortion as a free parameter does 
not improve the reduced \ 2 - The 2a upper limit on /i is 



fj, < 5.8 x 10" 



(8) 



This limit is an improv ement on the previ ous limit of 
9 x 10~ 5 set using FIRAS (|Fixsen et al.lll996h . Although 
it is numerica l ly sim ilar to the limit reported recently by 
iGervasi et all (|2008bf) . we believe it to be a more robust 
limit because we have allowed additional free parameters 
in our fit. 

5. DISCUSSION 

We have presented evidence for isotropic radio emission 
detected by ARCADE 2 beyond what can be explained 
by Galactic emission and the unresolved emission from 
the known population of discrete sources. The excess 
emission is consistent with a power law, with an index of 
—2.56, which is significantly flatter that what might be 
expected from an unidentified population of faint, diffuse, 
steep spectrum (index ~ —2.7) radio sources associated 
with star-forming galaxies. We have also examined and 
placed limits on two classes of spectral distortions to the 
CMB. Such distortions are not supported by the data 
and cannot explain the excess emission, as is illustrated 
graphically in Figure 4. 

It is possible to imagine that an unknown population 
of discrete sources exist below the flux limit of existing 



surveys. We argue earlier that these cannot be a simple 
extension of the source counts of star-forming galaxies. 
As a toy model, we consider a population of sources dis- 
tributed with a delta function in flux a factor of 10 fainte r 
than the 8.4 GHz survey limit of iFomalont et al.l (|2002ft . 
At a flux of 0.75 fiJy, it would take over 1100 such sources 
per square arcmin to produce the unexplained emission 
we see at 3.20 GHz, assuming a frequency index of —2.56. 
This source density is more than two orders of magnitude 
higher than expected from extrapolation to the same flux 
limit of the known source population. It is, however, 
only modestly greater than the surface density of ob- 
jects revealed in the faint est optical surveys, e .g., the 
Hubble Ultra Deep Field (|Beckwith et al.ll2006f) . The 
unexplained emission might result from an early popula- 
tion of non-thermal emission from low-luminosity AGN; 
such a source would evade the constraint implied by the 
far-IR measurements. 

We believe our result to be robust to several poten- 
tial sources of error. Underestimating the level of Galac- 
tic emissionjs a potential contaminant. As described 
bv lKogut et al.l (|2008[ ), however, the expected contribu- 
tion of the area around the North Galactic Pole is only 
~ 0.5 K at 1 GHz, with relatively tight errors. Further, 
the Galactic emission is estimated with two independent 
methods along several lines of sight. 

Correcting for instrumental systematic errors in mea- 
surements such as ARCADE 2 is always a primary con- 
cern. We emphasize that we detect unexplained emis- 
sion at 3 GHz with the ARCADE 2 data, but the result 
is also independently detected by a combination of low 
frequency data and FIRAS. The unexplained emission is 
detected above 3cr with any two of the three data subsets: 
1) FIRAS and low frequency radio data, 2) ARCADE 2 
and low frequency radio data, and 3) ARCADE 2 and FI- 
RAS. The result is therefore robust to problems in any 
one measurement. 

We conclude that the three most important sources of 
error, Galactic emission, instrumental systematic errors, 
and radio emission from the faint end of the distribution 
of known sources, are unlikely to be sufficient to explain 
the excess emission presented here. 

It is a pleasure to thank the staff of the Columbia Sci- 
entific Balloon Facility for their support of ARCADE 2. 
We thank the undergraduate students whose work helped 
make ARCADE 2 possible: Adam Bushmaker, Jane Cor- 
nett, Sarah Fixsen, Luke Lowe, and Alexander Rischard. 
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Fig. 4. — Fit of ARCADE 2 data, FIRAS data, and data from low frequency radio surveys. The upper plot shows (solid 
line) a fit with three components: a frequency independent CMB contribution, a power law amplitude, and a power law index. 
The lower plot shows the fit residuals. The dotted line shows the expected shape of a /i distortion. The amplitude of the 
plotted distortion is 100 times the 2a upper limit allowed by the fit. The dashed line shows the shape of a Yff distortion with 
an amplitude equal to the 2a upper limit. T he addition o f eithe r a fj, dis tortion or a, Ya dist ortion as a free parameter is not 
supported by the data. Data po ints are from lRoeer et aD (|1999h (cross). [Maeda et alJ 1 1999f l (asterisk), Has lam et al.l l|1981l ) 
(triangle), Reich & Reich] f|1986l ) (square), ARCADE 2 (diamonds), and FIRAS (heavy line), corrected for Galactic emission 
and an estimate of extragalactic radio sources, as shown in Table 1. 
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